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Abstract 

We establish a new framework of finite temperature field theory for Yang- 
Mills theories in the physical phase space eliminating all unphysical degrees 
of freedoms. Relating our method to the imaginary time formalism of James 
and Landshoff in temporal axial gauge, we calculate the two-loop pressure and 
provide a systematic and unique method to construct the additional vertices 
encountered in their approach. 
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INTRODUCTION 



This paper is devoted to establish finite temperature field theory on the physical phase 
space of nonabelian gauge theories. It is the intrinsic nature of every gauge theory that the 
whole configuration space contains gauge group orbits, and gauge transformations generate 
shifts along those orbits. Gauge equivalent field configurations are physically indistinguish- 
able, therefore only transitions between distinct gauge orbits contain physical information. 

In standard field theory the problem of superficial degrees of freedom is attacked by 
the introduction of a gauge fixing condition. However, due to Gribov's ambiguity JTJ], the 
gauge orbit space of Yang-Mills potentials cannot be parametrized uniquely by potentials 
satisfying a local gauge condition. A gauge condition surface in the entire configuration 
space contains gauge equivalent field configurations. 

Moreover, there are may exist certain field configurations where the gauge fixing surface 
is tangential to gauge orbits, corresponding to the zeros of the Faddeev-Popov determi- 
nant. Perturbatively, this entails, depending on the gauge chosen, unphysical poles in the 
propagator that have to be defined properly. In particular for the class of axial gauges, 
n ■ A = 0, = (1, n), the spurious poles at po = pn have to be treated with the co-called 
Leibbrandt-Mandelstam || prescription, which has to be modified H for lightlike. It 
should be mentioned that Landshoff's a-prescription M also gives the correct exponentia- 
tion in a Wilson-loop calculation up to order g 4 . 

At finite temperature in the real time formalism (RTF), analogously to the zero tempera- 
ture case, a temperature-dependent pole-prescription arises naturally within the framework 
of Hamiltonian quantization ||, and for the particular choice of temporal axial gauge (TAG), 
(n = 0), a RTF has been developed successfully ||. 

In the so-called imaginary time formalism (ITF), however, the energy can only take on the 
discrete Matsubara frequencies po = 2irinT '. Thus if one naively heats the unphysical degrees 
of freedom, an unresolved pole remains for the zero mode for momenta pn = 0. Luckily, in 
the particular case of the unresummed imaginary part of the transversal structure function of 
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the gluon self-energy, these factors cancel out for symmetry reasons and the straight-forward 
application of ITF Feynman rules works |7J . 

At the contrary, in TAG the naive ITF propagator contains unregularized singular factors 
1/po at zero Matsubara frequency. This problem has been circumvented in earlier works || 
by the ad hoc assumption that such poles have to be dropped. Although there exists no 
justification for doing so, the leading order self-energy is found to coincide with the results in 
other gauges. The deeper reason for this is, however, not the correctness of this prescription. 
In fact, in general axial gauge, the dependence on the gauge fixing vector n, and thus the 
axial poles, completely cancels out algebraically Jf] which gives the proof that no prescription 
enters at that loop level. In fact, a closer inspection of the corresponding expressions in TAG 
reveals that this cancellation takes place in that gauge too. Moreover, this must be the case 
since the leading-order expression for the self-energy is nothing but the hard thermal loop, 
and thus a physical, gauge independent quantity. 

It is well known, that the consistent calculation of the next-to-leading order contributions 
requires an appropriate resummation of propagators and vertices [[|. However, since the two- 
point function is a gauge independent quantity only on the mass-shell, one cannot expect 
the self-energy to be prescription independent off the physical dispersion relation. In the 
light of this line of arguments, the non-Debye screening behavior fllOUlll which contradicts 
results obtained in Coulomb gauge ||12|| , and using Polyakov-loop correlators |T^JT3|] appears 



to be rather an artifact of the off-shell calculation and not of the ad hoc pole prescription. 
Adopting the on-shell definition of the Debye-mass proposed by Rebhan |14] we expect that 



in a pragmatic calculation keeping the undefined l/po quantities those unphysical poles 
cancel out. 

Taking seriously the naive formulation, it was proposed quite recently [[HJ to regular- 



ize the divergent l/po expressions by a temperature dependent expression. However, this 
proposal may only serve to give the expressions an intermediate meaning, and it has to 
turn out irrelevant in the calculation of physical quantities. Moreover, it is not clear if this 
prescription leads a truly temporal propagator and can be adopted unambiguously without 



introducing ghost fields. Independent of any ad hoc method to get rid of the temporal pole, 
there is no justification for the naive application of ITF Feynman rules in TAG. 

It has been pointed out by James and Landshoff [[HJ already some time ago that the 
temporal pole is related to the free motion of the longitudinal modes of the gauge field. 
This violates periodic boundary conditions that are necessary to set up ITF. James and 
Landshoff invented a new formalism in which the longitudinal fields remain unheated and 
only the remaining physical degrees of freedom attain a temperature. Within this formalism, 
the longitudinal part of the propagator is automatically free from the 1/po singularity. It 
was argued that one can obtain the same answer for the two-loop pressure as in other gauges. 
The main drawback in that formulation is, however, that one has to construct physical states 
by explicitly solving the Gauss law which gives rise to additional time-independent vertices. 
This results in unwieldy expressions and it appears difficult to establish a resummation 
program. 

In the present paper, we advocate a different route to attack the problem. Based on a 
Hamiltonian formulation of the theory, we are able to eliminate all unphysical degrees of 
freedoms from the Hamiltonian by introducing an appropriate coordinate system in the space 
of field configurations which allows to make a unique distinction between gauge degrees of 



freedoms and physical ones |L7| • Within this approach one neither encounters any unphysical 
poles nor any explicit construction of physical states is necessary. It is straightforward to 
heat the physical degrees of freedom in the resulting non-local Hamiltonian. In order to 
illustrate the new method, we calculate the two-loop pressure and compare the result with 
the corresponding expressions found using the approach of James and Landshoff. Rewriting 
the physical fields in the basis used in their investigation, we give a general strategy to 
construct the corresponding Gauss law states to arbitrary order and calculate them explicitly 
to third order in the coupling constant. 
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I. CONSTRUCTION OF THE HAMILTONIAN IN THE PHYSICAL PHASE 

SPACE 



In this section we construct the physical Hamiltonian of pure QCD by eliminating all 
gauge degrees of freedom. The basic idea may be illustrated as follows. Consider a point 
particle moving in a plane and rotations around the origin as symmetry group. Then different 
trajectories are gauge equivalent if they can be mapped one to another by rotations. Of 
course, the dynamics appears simplest in polar coordinates, where only the radial coordinate 
has a physical meaning. The angular momentum generates gauge transformations and the 
corresponding canonical conjugate position coordinate, the angle, may be fixed arbitrarily. 

We start our investigation with the pure QCD Lagrangian 

L = -\j d 3 xTrF^F, u = -~ (F^F^) (1) 

where F^ v = d IM A u — d u A^ + ig [A^, A u ) , and the Yang-Mills fields are elements of the Lie 
algebra of the gauge groupQ. To go over to the Hamiltonian formalism, we have to determine 
the canonical momenta = 5L/5A' M = F 0/ \ The momentum conjugated to Aq vanishes, 
Eq ~ forming the primary constraint. The corresponding canonical Hamiltonian has the 
form H = 2 (A\ E l ) — L = (E\ E l ) + V{A L ) - 2 (A , Q) where 

V(A) = ±{F»,F») (2) 

appears as potential and Q = V l [A] E l = d % E l — ig[A l , E l ] is the so-called Gauss law, V l [A] 
being the covariant derivative in the adjoint representation. 

The primary constraint must be conserved during time evolution. This yields the sec- 
ondary constraints E = {E , H} = Q ~ 0, where we implicitly assumed the standard 
equal-time Poisson brackets \ A a (x\ E h (y) \ = 5 ab guiyS 3 (x — y). 

Since the algebra of the Gauss law closes, ^Q a (x), G b (y)} = g f abc 5 3 (x — y)Q c (x) and Q a = 
H} = -f abc A b o g c we conclude that there are no more constraints, and all constraints 



1 The hermitian generators of the gauge group are normalized according to TvT a T b 
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are of the first class. Following the standard Dirac quantization procedure ||19|| , the Poisson 
brackets may be deformed to eliminate the constraints E ~ at the operator level. Since 
the gauge orbits generated by E are shifts of A only, and leave the other phase space 
variables A\ untouched, we select the gauge equivalent configuration in the phase space 
which satisfies A = 0. This amounts to simply dropping the canonical pair A , E from the 
Hamiltonian. On this hyperplane, the remaining constraints Q become time independent 
since the Hamiltonian now commutes with the Gauss law. Those constraints generate time 
independent gauge transformations on the remaining phase space variables E % and A 1 , 

E n = QEQ-\ A n = QAQ- 1 - -(dQ)Q-\ 

9 

where Q is an element of the gauge group. It is convenient to formulate the quantized 
theory in a functional representation. Representing the canonical momenta by the standard 
functional differential operator E(x) — > —i5/5A(x) the Schrodinger equation takes the form 

H* n [A) = (- ^ + V[A]j * n [A] = E n *n[A] (3) 

and the wave function is subject to the constraints 

g^ n [A}=V[A]^ n [A]=0. (4) 



In the approach of James and Landshoff |16[ the wave functions are constructed explicitly by 
solving this constraint in order to be able to perform the thermal trace over physical states. 
Alternatively, one may eliminate the superficial degrees of freedom by reducing the number 
of field components. In analogy to the example given in the beginning of the section, we 
parametrize the field configurations in the unconstraint configuration space by an 'angle' u 
and the remaining coordinates A 1 in the following way 

A = UAU- 1 - -(dU)U- 1 . (5) 

Here U[u] = exp {iguj) is an element of the gauge group generated by the Lie- algebra valued 
angle uj and A % = e l a A a , a = 1, 2 are the remaining coordinates projected out by the operator 



e that is normalized such that e^e^ = 5 a p acts as unity in the corresponding subspace. The 
constraint (|j) transforms into 

s 

g u\ = U— V[A, cop' 1 = (6) 

which tells us that physical variables are independent of the angle u. The Gauss law itself 
turns out to be the canonical conjugate momentum to the 'position' variable to. 

In the potential part V[A] we may simply replace the original gauge field by its physical 
components, V[A] = V[A], since it is gauge invariant under the transformation flj|) . For the 
kinetic part in the Hamiltonian, we need an expression for metric entering implicitly in the 
definition of the inner product in (JJ). The lower metric components can be read off from 
the differentials, (The anihermitian covariant derivative is meant to be taken with respect 
to the fields A.) 



(5A, 5 A) = (5A\ 5 A) + (6u V f , 6 A) + UA\ V5u) + (5u, V f V5 



i-e. g a p = 4 4 = 93a = = 933 = Vitv *> and we defined the projected 

covariant derivative as V a = V*e^,. We shall also need the determinant of the metric, fi 2 = 
det(g) = det(P), T> = V^V — V^V, and the inverse components, which read g af3 = 5 af3 + 
V a P" 1 ^, g a3 = g 3a1 = -V at P -1 , fi- 33 = P" 1 . The kinetic term in the Hamiltonian 
reads in covariant form (5 / 5X l ,5 / 5X l ) = (^fi^d / 5X A g AB , /x 5/5X B ^ where we have put 
X A = (A a ,uj). When this operator acts on physical wave functions which by virtue of (||) 
do not depend on the angle u, the terms containing the momenta 5/5u vanish. Thus the 
physical Hamiltonian is obtained by simply dropping the 3-components in the metric and 
the Schrodinger equation (j3|) together with the constraint (|j) is equivalent to the reduced 
dynamical system described by the Hamiltonian 



(7) 



where the wave function as well as observables are functionals of the coordinates A only. 
Apart from the determinant \x appearing in the Hamiltonian, this expression was already 



found in |2(J using a different construction. The determinant resolves the apparently existing 
operator ordering problem which is due to the necessary inversion of the operator T> in the 
kinetic term. This problem was also discussed previously |2I|| , however without definite 
solution. The operator \i gives contributions ~ (5 3 (0)) 2 , which may be dropped if one is 
only interested in the local properties of the theory. 



II. FEYNMAN RULES IN THE PHYSICAL SUBSPACE 

The Hamiltonian derived in the previous section may serve to establish a set of Feynman 
rules in the physical subspace. Since it no more contains unphysical degrees of freedoms, 
it is straightforward to establish finite temperature field theory by heating the field A. We 
emphasize that no explicit choice of the projection operator e l a corresponding to particular 
coordinates is necessary so far. However, (|7|) contains a non-local operator which makes 
the theory unwieldy to deal with. Alternatively, one may introduce an auxiliary field and 
rewrite the Hamiltonian density in the following manner 



H = }-E a ' a1 E a ' a + - 



§ a i(y a E a ) a + (v a E a ) at $ e 



2 4 r H ' 

(8) 



The electric field E is the canonical conjugate to A in operator representation, 
[E(x) a ' a ,A(y) f3 ' b ] = -iS ab S af3 S 3 (x - y). The last term in <E) just subtracts off the 

xo=yo 

trace of the operator T>, which amounts to drop all Feynman graphs which do not contain 
at least one $ a (V a E a ) a vertex. 

For the particular choice of purely transversal fields, d % A 1 = 0, the operator sandwiched 
between the covariant derivative and the electric field becomes the spatial transversal pro- 
jection operator, e l a e J a = 5^ — d l d^ /d 2 = T u (<9). In that case, the $>VE terms in (H) turn 
into a single three-vertex <&AE and not two-vertex remains. We observe that even for 
a non-transversal choice of A, E, the perturbative E L E L propagator is compensated by the 
E L § two-vertex and the Green function. This corresponds to the fact that in pure QCD 
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zeroth order longitudinal states are pure gauge degrees of freedoms. We want to point out 
that although all fields are purely transversal they do not coincide with the transversal com- 
ponents of the original gauge potentials, and the transversal choice must not be confused 
with the Coulomb gauge. 

Splitting the Hamiltonian into a free part 

% — - ^E a,i T ij E a,i - A i ' a /\T ij A i ' a + $ a A$ a ) 

and an interacting one 

u u$AE i u&A i uV(A) 

itww — it t rt + rt v , 

r H® A = n f abc <& a A Ti ' c d i <& b — —f abc f ade (^ b j\ Ti ' c j\ T i' e ——^) d 

n 2 

, r ( V ( A ) — _gf abc Q3 J)Tt,a j^Tj,b j^Tt,c _|_ y_ jabc jade j^Ti,b j^Tj,c j^Tt,d j^Tj,e 

the former one gives rise to the propagators 

A Ti A T A(p ,p)=-T i i(p) 1 



9 -*0 

Po~P 



2 



o -»9 

Po~P 2 



E^E^)(p ,p)=-T^(p)-/- 
= ($$)( po ^) =-4- 

and from the second one we shall only need the expression for the $ a A Tt ' b E T i' c vertex given 
by (-l)gf ahc 5 ij . 

For two reasons, we do not express the electric fields in terms of time derivatives of 
the gauge potentials. Firstly, this would again involve inverting non-local operators, and 
subsequently render the theory untractable. Secondly, the kinetic term and the potential 
term are both by themselves physical observables. Keeping the E fields, one preserves the 
possibility to calculate the electric and magnetic dispersion independently. 

At the present stage, the propagator of the field $ only serves to write the inverse 
Laplacian in a way convenient for calculating quantities in ITF and remains unheated. In 



ol = 



o2 = 



a3 = 



FIG. 1. Graphs contributing to the two- loop pressure in the physical subspace. Full lines 
correspond to A propagators, the wavy line to the E 1 - field propagator and the broken line to the 
auxiliary field. 

the due course of an eventual resummation, however, consistency may require to assign the 
auxiliary field a temperature dependent Green function. This does, of course, not contradict 
the original nature of an auxiliary field, since to a given order, the perturbative inversion of 
T> in the non-local Hamiltonian is not unique. 



III. THE TWO-LOOP PRESSURE 



wc 



In order to compare our method with the construction of James and Landshoff [|l| 
calculate the two- loop pressure given by the diagrams depicted in Figl. We observe that 
apart from the third diagram, only transversal fields contribute to the pressure. 

The temperature dependent contribution reads {Cn = N(N 2 — 1) for SU(N) ) 

where Z = (j\ + cr 2 + cr 3 and <7j corresponds to the three graphs in Figl respectively, 

1 (3-z 2 ) . . 1 z (l + z 2 ) . 

° 1 = l pq ( npUq + Hp) ' ° 2 = 2 (p + q) 2 ^ pnq + Hp) ' 

1 1 (1 + z 2 ) . 2 2 , . . 

z = pq/pq, p = \p\, n p = (exp(p/T) — l) -1 , which gives the correct answer 

p(2) = L Z W = 9 2 C N T A 
V 144 ' 
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FIG. 2. Graphs contributing to the two- loop pressure in the approach of James and Landshoff. 
The blob denotes an additional three-vertex due to the explicit construction of Gaus law states. 

Note the simplicity of the calculation. As opposed to that, the approach of James and 
Landshoff involves zero-temperature longitudinal propagators and heated transversal ones. 
The calculation of the pressure calls for the calculation of 12 graphs depicted in Fig2. 

In addition to the usual three- and four- vertices of QCD, a static vertex denoted as a 
blob in Fig2 enters from the explicit construction of physical states satisfying the Gauss law 
constraint (Q). The (TTL) part of that vertex is asymmetric in the transversal legs which 
we indicated by an asterisk. We neither want to repeat the details of the formalism nor the 
lengthy but straightforward calculation but rather state the result. Using the notation of 
Eq. (H) we find Z = J2l=o K i where (The contributions n\ . . . k 4 have already been calculated 
in PI.) 

„2 , 2\ 



1 1_ (l + z z 
' ' 2 pq (p + q) 



{n p n q + n p )(p 2 + q 2 ) + (2n p + 1) ^ P 9 

p 2 + q l 



K5 ~ 8T (p + q)^ Up+ 2 h 

K = 02, K2 = ^3 = «2, K 4 = Cr 3 , K 7 = K 6 = K 5} Kg = ~K 5 , Kg = 0, K 10 = G\. All 
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(TLL) graphs cancel. Furthermore, since K-y + K2 + K3 = only the term K4 remains in the 
(TTL) contributions that coincides with the third graph in Figl in our calculation. It is 
interesting to observe that only the (TTL) part of James and Landshoff 's new vertex plays 
a role and that the only non-vanishing contribution comes from diagrams which contain a 
pair of that vertex part. 

In the case of the two-loop pressure, one may even replace the longitudinal .E-fields in the 
original Hamiltonian by the vertex part of the Gauss law constraint, i.e. dE — > gf abc E b '' i A c '' 1 , 
and drop the longitudinal A fields to get the correct result. However, it is clear from our 
investigation that this oversimplifying guess does not give the correct answer in general. 

IV. CONSTRUCTING GAUSS LAW STATES TO ARBITRARY ORDER 

There are two different ways to handle gauge degrees of freedom at finite temperature. 
Based on BRS invariance, Hata and Kugo |23| constructed a theory where the Boltzman 
factor in the thermal average gets replaced by exp (mN c — H/T), N c being the ghost number 
operator, and the trace is expanded to include ghost fields and all degrees of freedom of the 
gauge potential. They demonstrated that with this weight thermal averages of operators 
corresponding to observables are the same as in the projected ensemble. The advantage of 
their approach, which is the standard way finite temperature field theory is handled, is that 
all degrees of freedom are heated which results in simple Feynman rules. However, in the 
particular case of TAG, this construction does not work || and one has to go back to the 
projected ensemble involving physical states only. 

There, the thermal average of an observable Q is defined by 

(g) = Z- 1 ^(Phys|e-^ T g|Phys) 

where physical states satisfy the Gauss law, Q |Phys) = 0. In the formalism of James and 
Landshoff those are constructed by acting with a unitary operator on the free transversal 
states |T), 
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|Phys) = i?|T>, R=Y,9 n Rn (10) 

n=0 

where R n was determined by acting n + 1 times with the Gauss law operator on the physical 
states. There the longitudinal components of the gauge potential correspond to gauge de- 
grees of freedom. In our approach, the decomposition ([5|) allows us to identify the longitudi- 
nal components with the 'angle' variable, A L = du, whereas the remaining field components 
read 

A T = A — A L = e l9LU Ae- l9UJ - - (de i9 ") e~ i9UJ - du. (11) 

9 v ' 

Recalling that the field A contains but physical degrees of freedom, one realizes that the 
transversal components of A are physical only to zeroth order in g. Conversely, the choice 
to keep the longitudinal components A L unheated to all orders corresponds to the fact, that 
the physical field A is transversal only to lowest order, which in turn means that there do 
appear higher order heated longitudinal modes contained in the states that satisfy the Gauss 
law. 

The operator A can be expressed in (A T , A L ) coordinates by virtue of (|TTJ) 

00 / nn i \ 00 (ia) n ( n ) m n 

A = Y. 9 n A n = A + {du) + -d) = £ M- [ u, A + ^-du] (12) 

V 9 J n\ n + 1 

(n) (n-1) (0) 

where [ X,Y] = [X, \ X,Y}}, [ X,Y] =Y denotes the multiple commutator. A is the 
counterpart of the Gauss law operator to arbitrary order, where the first few terms read 

A = A T + ig[u, A T + ±8uj] - [u, A T + ^duj]] + .... 

Since the change of the basis from physical to transversal states is mediated by a unitary 
transformation, the corresponding operators A and A T are unitary equivalent according to 
AR = RA T . Collecting terms by orders of g, this leads to the recursion relation 

n-l 

[Ao, R n ] + ^ An-mRm = 0. 
m=0 

The strategy to solve this recursion may be motivated by the following observation. R n is 
given by the action of an (unknown) operator on the sum which inverts the commutator with 
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A = A T . Recalling the equal time commutator [A Tl ' a (x), E T ^ b {y)] = i8 ab T l i(d)5(x — y), 
one may, roughly speaking, construct R n by 'multiplying' the sum with E T . 

In particular, for n — 1 one has to study the equation [A Q , Ri] + A\ = that has the 
solution 

R x = i(E T ■ ii) := i J £zE TA (z)A\(z). 
Plugging in the expression for A\, the explicit form of R\ can be written as (ou = dA L /d 2 ) 

Ri = it bc j d 3 z [^dA L - a y z ) E T ' b (JF^z) + \A L ' c {z)^ 

which coincides with the result of James and Landshoff. 
The recursion for n = 2 reads 

[Ao,R 2 ]+A 1 R 1 +A 2 = 0. (13) 

Guided by what we have learned above, one would naively guess R^ 1 ^ = i (E T ■ A 2 )+ 1 / 2-R1-R1 
which when commutated with A gives [Aq, R 9 2 uess ] = —A 2 — x / 2 AxR\ — l / 2 R\A\ that cancels 
the third but not the second term in the recursion ( |i~3"D since A\ does not commute with R\. 
We therefore add a suitable chosen term proportional to the commutator which compensates 
the wrong order in the R\A\ contribution. The solution of (|i~3|) reads 

R 2 = i{E T - A 2 ) + ^R 1 R 1 - 1 -{E T ■ [R u At]) 

which is unique up to operators that commute with A . 

One may continue further and calculate the n = 3 contribution to R, 

R 3 = ^Rl + ~ (V- [ R 1} AiYj - % -R x (E T ■ [R l7 A,}) + iR x (E T ■ A 2 ) - 

-1 (E T -[R 1 ,A 2 })+t(E T -A 3 ) J 

and a calculation of higher order terms proceeds analog similar lines. We note that our 
construction has formal similarity with the Foldy-Wouthuysen transformation in quantum- 
mechanics. 
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Our result does not match the R2 contributions found in [22] that also contain time 
derivatives of the longitudinal fields, but agrees with the argument given by James and 
Landshoff that those time derivatives should not be present in R. We also find that the 
exponentiation conjectured in cannot be confirmed. Although the Ri terms appear 
with the correct factors, a complete exponentiation is spoiled by an increasing number of 
commutator terms, which is consistent with the nonlocal Hamiltonian (0) that also contains 
an infinite number of vertices. Unlike as in pure QED, where the radiation gauge eliminates 
all gauge freedoms from the Hamiltonian in a local manner, and where R does exponentiate, 



is has been argued [20| that in nonabelian gauge theory no such canonical gauge exists. 



V. CONCLUSION AND OUTLOOK 

We formulated a finite temperature framework for pure QCD which is based on the 
elimination of all gauge degrees of freedom. In contrast to the former approach of James and 
Landshoff, who explicitly constructed physical states by solving the Gauss law constraint, 
we eliminate spurious degrees of freedom at the operator level which allows to heat the 
remaining degrees of freedom in a straight-forward way. We do not encounter any pole 
ambiguities which exist in the naive imaginary time formulation of axial gauges. 

Although our effective Hamiltonian contains a kinetic term non-local in the fields, it is 
possible to find a local formulation by introducing an auxiliary field. We compared our 
theory with the construction of James and Landshoff for the particular case of the two- 
loop pressure and found that the number of Feynman graphs is reduced drastically in our 
framework. Since our construction allows to make a clear distinction between physical and 
gauge degrees of freedom, we can also line out a strategy of how to explicitly construct Gauss 
law states to all orders. Those are calculated explicitly up to third order in the basis of the 
free transversal states of the gauge field. 

Clearly, we will establish a resummation program in our formalism. Since the Hamilto- 
nian only contains physical observables, the location of the poles in the propagators contains 
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intrinsic physical information on the the on-shell dispersion relation. It would be interesting 
to compare with results obtained in the usual approach, where gauge-fixing independence 
of the poles has to be and was proven . 

Furthermore, we only dealt with pure QCD, which by construction excludes the calcula- 
tion of the Debye-mass that would require the gauge-invariant inclusion of a charge density 
in the Gauss law constraint. Conversely, if it turned out to be possible to include charges 
in the present formalism it would be possible to separate effects from external charges and 
those induced by pure QCD. We are going to investigate on this question. 
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